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■ METHOD FOR BISPECTRAI, PICKING OF ANELLIPTTrar. rjMn ' 
•• ; CORRECTION PARAMETERS . „...,:- 

The field of this invention is seismic prospecting 
The invention relates particularly to processing of 
seismic traces in a common midpoint gather. 

The invention. is more precisely related to a method 
5 Pf . determining velocity V and. anellipti city ^parameters 
... necessary to perform processing including NMO correction 
of seismic traces. . , ■ ' ... < :-: 

Seismic prospecting usually, consists'," oi emitting 
K seismic waves in the ;sub.s6il . using. 6ne or' several seismic 
10 sources, and making surface records- of seismic data 
corresponding to seismic; waves, reflected on geological 
.. interfaces in the subsoil as a function of time (also 
, called, reflectors): using receivers (also called geophones 
or hydrophones) depending on,. whether the prospecting, is 
15, being done on land or. at . sea) and then' processing these 
data to extract useful information about the geology of 
the subsoil. 

The record of the seismic energy produced by each 
receiver during data acquisition is called a seismic 
20 trace. 
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One conventional • technique for seismic prospecting 
is multiple coverage, in which sources and receivers are 
arranged such that several seismic traces are grouped 
together at the same midpoint (in other words the point 
at equal distance between the source and the receiver for 

a given trace) . 

Although seismic traces contain useful information 
about seismic reflections and geology of the subsoil, 
they also contain noise components. 

One of- the first-/ purposes) of the processing of 
seismic data is ta eiirhihate/.- o^at least - attenuate, 
these unwanted noise components such that the useful 
information can be clearly identified and interpreted. 

One conventional method used to attenuate these 
15 noise components is Common Midpoint (CMP) gather. Traces 
with the same midpoint are then grouped as a function of. 
the. distance between the source and the receiver (called 

, the offset). " • ' "'■'■f%%. 7 : '. . '". . 

,.; :■))''■' i n " general , ; ^representation 'of seismic dataVin4 images;: 

2 0 requires the ■ use. [ of 'proces sing : including : 7 

- - a so-called I TZO - . (Transform • to. Zero Offset) 
operation designed to compensate for the NMO effect of 
pathsj by ■ correcting the ; arrival^ time ! of reflections: to 
. the arrival time of traces with, zero . off set ,| v-, 
25 . - and a . .'. -. migration 7 operation- designed : tpl reproduce 
..the. correct shapes, of, the geological interfaces'. -/ ; : - / ; . • 
. Although these TZO and migration operations are 
usually carried out one after the- other, they can be 
carried out jointly. This is the case particularly when a 
time migration is done before stack (PSTM - Pre-Stack 
Time Migration) . 
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in a simplified manner, the TZO operation simulates 
acquisition of seismic data by sources and receivers 
placed at the common mid-point. 

The purpose is to., add records illuminating the same 
1 ' P ° int in - the subs °il • to increase, the signal to ; noise 
ratio-, and ,, the ., primary reflections. .. to secondary 
reflections, ratio, and thus benefit from "the advantages 
of "multiple coverage"., 

in order to make a zero offset image, a so-called 
obliqueness correction method or NMO . (Normal . Move. Out) 
, s correction method is used. . ' ""'*'*•. :,,< 

If it is. assumed that the subsoil is 'horizontally 
stratified with no lateral variation of the propagation 
velocities, it can be shown that records that have the 
property of illuminating the same point in the subsoil 
are records that have the same midpoint. , 

, . However,, the image of a reflection in the ' subsoil 
• ai " riVeS at variabl 1, times depending : on the offset. 

Therefore,, in ?^|p>tack reflections, it is necessary f 
, t0 - art : >^ cor^ebyng the ; various} records \t0 %^the m 
' all to a common zero offset reference; ^} • Z ^^MS^tf' 1 :' V 
; Historically,,: ^j^^^^f^j^^^^^^ 
particularly simple model of the subsoil; a homogeneous 
model with horizontal reflectors . 
,., In : > this r< ^el ? ,. reflections;' associated with i each 
subsoil reflector; are. theoretically ! aligned ^ along 
hyperbolas also called; indicators, \ centred vertically 
above the midpoint." 

The arrival time of a reflection is then a 
hyperbolic function of the source -receiver offset, the 
shortest time being obtained at zero offset. 
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in order to. stack records in each gather, the NMO 
correction straightens the hyperbolas until they are 
theoretically horizontal. 

The NMO correction is then made based on the 
5 following hyperbolic equation for the travel time t after 
reflection, associated V with . a; source-receiver offset x 
pair: ■ 

in which to represents the zero offset travel Vime; and V 
10 denotes the average propagation velocity of waves in the 
subsoil . 

'The simplistic model mentioned above is based 
particularly on assumptions such as low angles of 
incidence and an isotropic medium. 
15 But assumptions made for' the simplistic model are 

• too restrictive: to ' describe a . complex medium, and in 
, • particular cannot ' be "applied ..to v ' the.' propagation of ; 
seismic waves in/ an anisotropic medium (medium iri;. ; which 
'the' velocity of ~ ' waves ' can " vary" ," depending on the 

20 propagation direction) . ^ ..\ ( ,. /. 

; ; ' The' use-'of'-'W less simplistic mode? 'is necessary, 

.• particularly due to: . "{." : \ '• v ,.,/../': ( 

- the use of long streamers, for acquisition in deep 
offshore, which in particular means that- traces "with far 

25 offsets are acquired; 

_ observation of ' anisotropy " in ." clayey type 

sediments . 

The generally accepted assumption consists of 

modelling an anisotropic medium as a stack of 

30 transversely isotropic layers with a vertical axis of 
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symmetry. This is known as Vertical Transverse Isotropy 
(VTI) . 

Thus, it has been proposed to determine NMO 
corrections to be made by: 
5 - introducing vertical non-homogeneity into a 

homogenous VTI medium model, as presented in the document 
"ALKHALIFAH T. arid .TSVANKIN I. , 1995;. Velocity analysis 

for transversely isotropic media:- Geophysics^ 60, 1550- 
' 1566"; ■ • ■ .,, • 

10 . or by •; '• ■ .••v-v; '<-\ .. 

- introducing VTI anisotropy into a model with 
) stratified isotropic layers, as demonstrated in document 
"SILIQI R. and BOUSQUIE N. , 2000, Anelliptic time 
processing based on an offset hyperbola approach, 70 th 
15 Ann. Internat. Mtg. : Soc. Of Expl . Geophys., 2245-2248". 

This second approach combining vertical non- 
homogeneity and VTI anisotropy to give a new model of the 
subsoil, appears to be the best approach in most real 
. .[.; .cases studied.;.;.. - * ' y : " 

The following .<^j^oh^ 

't after reflection, ' using ; an anelliptic "offset 
hyperbola, associated with a source - receiver offset x 
pair, is derived from this model : ' - ' "' V 

t (V,7l) - 'J±Jf . \[ to V X 2 , .. .. 

, : v " ; W 0+ AT^) + ^-^ Ration (la, 

" where V is the - .velocity: conventionally used in seismic 
corresponding to small offsets and V is, a- parameter, 
called the anellipticity parameter. 

The document "SUAUDEAU E. and SILIQI R. , ,2001, 
Anelliptic pre-stack time migration, Annual International 
Meeting, CSEG Expanded Abstracts" also proposes to 
include the NMO correction by anelliptic offset hyperbola 
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in the equation for the path time used during the PSTM 

migration operation. 

The PSTM migration equation is conventionally 

expressed in the form of a double anelliptic offset 

hyperbola >. which is the sum of two square roots (the 

Double Square Root (DSQR) equation) . : . " 

,.0- Taking account of anellipticity, the expression for 
this constant offset equation, becomes : ."• ' ■ 



t = 




(l+S^V 2 VU+8tiJ (l+8n)V 2 



10 Equation (lb) 

where : 

- the V and t) parameters are those mentioned above, 

- Xu is the coordinates of the midpoints, 

- x - Xn, represents the migration aperture, 
15 - h is the half source-receiver offset, 

- t 0 is the 'double time , for zero aperture of the 

■ : operator . ;' \ / . ; ■ .-. .- : ■:, .-. v : }/s^II\s:\ ;a' .^.iL-'-'v •■ •>.•.. 

}.} Note, that when the aperture x - ^ of the migration 
" is zero'/ - the PS^Vcorrection^ equation ^ becomes the 

' 20 NMO correction; 'ea^ation ^ . .'Thus, the NMO. correction 
forms a special case of the PSTM migration; the case of 
PSTM migration with . zero aperture . 

Therefore, finally;, in order, to perform seismic data 
processing including NMO correction taking account of 
25 vertical heterogeneity and VTI type anisotropy, it is 
necessary to determine the velocity V and anellipticity r\ 
parameters . 

• The said parameters V and n can conventionally be 
estimated by carrying out two passes such that : 
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- during the first pass, the distribution of 
velocities V along the time "axis is estimated, only using 
data with near offsets; 

- during the second . pass, the anellipticity r\ is 
5 estimated along, the time axis, using: : - • 

- the distribution ; : of velocities : . determined 
during the first pass, and;/, ., v .^/;., 

- all data (including; data j with^f arf of f sets) ,;-;;h 
\ The. document entitled: >\ ;tt:> : : : u Z 

10 .' •■ "SILIQI . R . •/ . . 2p01, \;>Technolog;ical : * leap 4- in - time 
processing focuses the data throughout, anisotropic media:. 
■■. First Break, .19,. No. 11, 612-618 > :yv 
has also demonstrated that the parameters V and *n carl be 
estimated in a single pass, during which bispectral 
15 analyses are carried out in order to simultaneously pick 
the V and T| parameters along the-;, time axis, using all 
• data. : "■ ' . ■ : .V.-: . :; . ' " -v ;..-.;fyA v •" -\ . ' ' a-: : " .. ' ; 

However, dense analyses of NMO correction parameters 
are' preferably:.- made?- wheift the correction": \ equation '"no ■ 
20 longer depends on the., time t 0 (we" then talk about static :; 
: • NMO correction) '. W\^st . - '^0 f $r$ /W^f' : dl : . '$' : ' : ?% : 'W : ilT- " - f0- 
A static correction effectively provides a- means- of 
offsetting all samples- each forming; traces; ; of the , same : 
time St,; for a given, of f set . 1 : ' * \ -^'^ — . ; 
25 Thus , \ when; ' a ' static : correction' is made, 'the "number 

of calculations - to be made ' can s t>e ; '' significantly • reduced 
and the trace stretching phenomenon is eliminated, which 
makes the said dense analyses viable. 

At the moment, the only way to make velocity and 
30 anellipticity analyses is to use two passes, in 
particular making use of parabolic approximations of NMO 
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residuals, in order to obtain a dense picking of the V 

and r| parameters. 

In this context: 

- the first step is to estimate residual velocities, 
5 using near offset data, according to. a first estimate of 

velocities " v " ;■■ ;"• * A : :t *V. ;■ •;■ 

- the next step is to estimate anellipticity on all 
data, using velocity ^updates dpne^ previously, o 

A mute law must also be defined to. estimate residual 
10 velocities,^ so as: to only keep data that are considered 
■ as being . near offset data* '^v V ; ■ v. vc 

However, the efficiency of the two-pass analysis is 
particularly, sensitive to the choice of such a mute, law. 

On the other hand, the anellipticity T| is estimated 
15 mainly making use of far offset data. 

• But, the anellipticity est ima.te made for far offset 
data is not very precise, . such that the correction made 
• is .finally: inaccurate : - G&- V. ^ ■ / V\- -v-.^; rV,> ; 

One purpose; pf,v.the$ invention, is • to ^ overcome these 
20 limitations and disadvantages by proposing a processing 
including a static NMO correction which is more efficient 
and more " precise „vV - v jif^i * - . \ . , U . 

More, precisely, the purpose of the invention is a 

"dense determination of velocity; V and _ anellipticity . t\ 
25 •;. parameters in a single pass; using- all; available data/ in 
: other words; making- use 'of all offset ranges 

To achieve^ this , : the invention proposes a method of 
determining the velocity V and anellipticity r| parameters 
for processing seismic traces from a common midpoint 
3 0 (CMP) gather including an anelliptic NMO correction, 
comprising: 
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- a preliminary step to define a plurality of nodes 
(dtn, t 0 ) , the said nodes being indicative of parameters 
dtn and t 0 representing the NMO correction for the 
maximum offset and the zero offset travel time in 

5 hyperbolic coordinates, the. said preliminary step being 
followed by. . . Y . . , \ tl ■[■■'■ 

- for each node (dtn, * x 0 ) defined in the preliminary 
step, steps for the NMO correction, of traces in the CMP 
gather as a function of the values , of the said parameters 

10 dtn and x 0 . at the. node considered, and for calculating 
the ; semblance , function associated , with . the said NMO 
correction for the node considered; and... 

- for each picked time t 0 , a step including 
determination of the maximum semblance node (dtn(t 0 ) , 

15 to (t 0 )), 

- and a final step to convert the dtn(t 0 ) and t 0 (t 0 ) 
parameters so as to obtain the velocity V(t 0 ) and 

. anellepticity r| (t 0 ) laws.^V*i"v : C> ^.^ , . - . ^ ..v-: v : ; .• 

• Cv, According ; to a , first ^embodiment i of •■ the/ invention, . 
20 the processing "done, is a static NMO correction of. seismic 
' traces. '' ; ■ l ! S\ U S< , 

. According to a second . embodiment of the . invention, 
the^ processing done i is.; a . PSTM ■ migration of- seismic 
traces , the said PSTM migration . including a static . PSTM : 
25 NMO correction of the said seismic traces. 

' A preferred but non- limitative .aspect of the process 
according to the invention relates to the definition of 
the dtn and x 0 parameters relative to the velocity V and 
the anellipticity r| so as to make static NMO corrections, 
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according to dtn = — — 1 0 + 
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+ ■ ™* and T 0 = 
(l + 8r])V 2 



t 



o 



1 + 8ti 

Other aspects) purposes and advantages of this 
invention will become clearer after reading the following 
5 detailed description with reference to' the attached 
Figures in which: - ■ ; 

- Figure la shows the anelliptic offset hyperbola 
used to make the NMO correction and illustrates '■ the 
meaning of the parameters x 0 and dtn; ' - ! . ' v. 

10 - Figure lb represents the DSQR equation of the 

double anelliptic offset hyperbola of the PSTM migration 
and illustrates the meaning of the x 0 and dtn parameters; 

- Figure 2 illustrates the effect of the t 0 
parameter of the reflection curve corrected by the NMO 

15 correction;. — -■ [ - >s- : "=?s-.V 

■ v - Figures 3a . and 3b represent the analysis, volume 
(t 0 , dtn, t 0 ) in which the bispectral picking of the dtn 
• and x 0 f parameters' according to the invention is made ; 

•' . . - Figure 4 combines each of : (V, , V^). and (dtn/ . x 0 ) 
2 0 approaches, representing their corresponding bispectral 
analysis panel; . " 

- Figure . : 5 .. shows .; the r correspondence- between the 
(dtn/ t 0 ) and (V, Van)- pairs of NMO correction parameters; 

- Figure 6 shows a CMP gather of true seismic traces 
25 before the NMO correction and the bispectral picking of 

the (dtn, - t 0 ) parameters corresponding to this gather of 
traces; 



10 



- Figure 7 shows the semblance function and velocity 
and anellipticity functions deduced from the picking of 
the dtn and t 0 parameters in Figure 6.. 

- Figure 8a is "a flowchart representing the steps' of 
.; a : first' "particular 'e^pdiment: [ of ^' the' inyentibri;'' namely 

determination ' of" the Tv ancT r\ : 'parameters ' to * make an ".'nmo ; ". 
•correction; " ^:Vo V^ J 

*y '-"Figure 8b is a flowchart repre^enVing'the steps of v 

a second particular embodiment 'of the 1 invention, namely 
. determination- of -i** V and r\ parameters to- do a PSTM 
' migration ; ' ■ ; ' : ' : : " : ■ .« .. '■*. '= -%* ; ■ v^|£^v^;' ; '■ ^ ' ; I • --.r^r:': \ 

- Figure 9 is a diagram illustrating the different 
operations carried out in order to determine the V and r| 
parameters to make a PSTM migration. 

In general, the method according to the invention is 
a method for processing records' of variable offset 
seismic traces, and'; this processing' uses the recorded 
seismic traces to build common /midpoint (CMP) trace 
gathers , and the traces" in each: gather 'are subjected to 
an NMO correction. , « ^ — - i - , > 

•"' • "'' • : 'y-ii . • ; •.• vv\^'.:':- : v : .-''- : • ... s - ,v : -;- - : ; . V'.'-' "■ rf : s '- v f ' ;, ' : '-> : - •.• • 

; In particular, ther method according to, the invent ion 

1 ' i - •" P c '' .... .? •' ' '' • ^ ■ • >■■ U v a'>: ;■: - , v , 1 ••■'.■«-*., ■• ? - ■'■ . . - . t-'.*'; • >.••■« ■ 

.• •• *.*■•- ,\'V..V: •.'•*.••'>..*. V'V:.-' -'Vi ! - % 'v':^:-': v&r?'/..V^ s v^ •■• :,i >^ •'••'•"■<••..••■••?.''.•■•'>,.■• • •:• 

determines, the velocity y,. and /anellipticity r\ parameters, 
to perform processing including such an NMQ correction of 
the seismic traces of a CMP gather. 

For example, the said; processing may be : 

• • . , v7 ■■ ■: , - ;^v^v. ; >-, " -:V".--^- . . ; ".. - r 

: - an anelliptic NMO correction of the seismic 

traces'; . ' V - ' . .. 

- an anelliptic PSTM migration which, as we have 
seen above, jointly applies TZO and migration operations 
(in the following we will refer to a PSTM NMO 
correction) . 
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The following description is more specifically 
applicable to the NMO correction. However, considering in 
particular the two particular embodiments of the 
invention that will . be; described, below, we can see that 
this description is. equally. : applicable to. any processing 
including an NMO .: correction, particularly.. . processing 
including a PSTM NMO correction. " ' l -^J- , ^ 

. Two new .. parameters :.- dtn, and x 0 are, considered in 
order to, specifically determine y. and rp ■ >• ,c . 
:v xi ~ T o.;tha.t ; represents. \the. zero offset travels, time h in 
"hyperbolic coordinates" ,( see Figure 1) , ' ■. >• 

To = — — 1| : Equation (2) 
l + 8t| 

and, 

- dtn that represents the NMO correction for the 
largest off set. Xmax (see Figure .1) , ,. 

, , ,dtn = t x=x ^ . - t xM) namely: ... . fc , ... 



dtn ^ ,-^'V 0 +. 





4c " 2 










0 note: tl 


aab dtn 



— Equation (3) 



to* th^v.velocity;\V?Y.a^ :- the-:;aheliip1:icityVTi /while * x 0 is ; a 
perfectly" anellipticai parameter defined relative* "to the 
anellipticity r\f independent of V.; 



Accordingly, the 'velocity ,V . .and anellipticity r| 
parameters may. be calculated, in accordance with equations 
(2) and (3) , using the; following conversion equations (4) 
and' (5) : ■ ' ' '"' 



V = Xma Equation (4) 



12 



10 



and r\ = -I — -1| Equation (5) 

8 I 7 n 



Figure la shows, the . anellipt^ .hyperbola 
used to make the - NMO . correction and illustrates ' the 
meaning of the T 0 - and dt^ 

^ ? • The .said " hyperbolic^ r coordinates " are shown, in this 
Figure 1. Their origin on the time axis is taken at the 
intersection of,, the. said time, axis with the. asymptote 
tangent to the said far offset offset hyperbola (see. 
equation (la)). 

Using the (dtn, x 0 ) parameters,, equation (la) for 
the anelliptical offset hyperbola becomes: 



, 2 dtn(dtn + 2T Q ) 2 « ^. , v 
t = to.- x 0 +, /Tq+ — — x Equation (6a) 



k max 



15. . ; .The (dtn, x 0 ) parameters , defined f or the- yeloci t y V 
; ^ancU. anelliptici.ty , , t| •; area, thus v; used • ' to-, make .j;. the NM0:> : 
•correction CORR^- = t'*- to be applied ^ to * of f set x - 
:• traces " :; iridependeht;.:-of^- to V^"*^ : i>v . "-v^-K' -vi 

20 Therefore, this is a static NMO P correction. In other . 

.words, the data recorded on a given offset trace will all 
. be : corrected in the ; same , way ^f or, a^ 

independently of the time at which these 1 data were 
acquired. 

25 Consequently, the estimate of the velocity and 

anellipticity parameters is not disturbed by stretching 
of traces, usually observed when dynamic NMO corrections 
are made. 
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Similarly using " the two parameters (dtn, x 0 ) , 
equation (lb) for the DSQR double offset hyperbola for 
the anelliptic PSTM migration becomes: 



t = t 0 - *0+ 



|T; Ydtn(dtn4-2T 0 )(x-x m +by 'V^ 

1 4 i'%. : 



w.v-:-.--v., ■« .-i... fev+ ^ggT^g 2 < ta^^^<: 7 ? figyg^pn (. 6i >) ; 

'■' V V>;. ' *W? .» -," . A 

Figure lb shows the DSQR double offset hyperbola of 
the PSTM migration defined with, the (dtn, x 0 ) parameters. 
10 In the context of the PSTM migration, represents 

the maximum offset and aperture of the migration. . 

The (dtn, x 0 ) parameters defined relative to the 
velocity V and anellipticity r\ are thus used to make the 
PSTM NMO correction CORRpstm - t - x 0 to be applied to the 

15 offset x traces independent of x 0 (equation 7b)): : 

- ", . ' " . 1 . , "iu,'^-? '.i.-' t - ' : - * • 

Corrpsra (X). = -Tp + . K ■ dta(din-t2T t Xx-x+Mf; + & , dtn{dtn + 2f „)(x - x - hy~ , 

•": WV '' t VT : — , s VT" ~ ' • ; . ; ' 

Therefore, . this., is a static PSTM correction.'; In 
other words, the se ' of ' aperture samples x -"•'x m -'of'' an 
"iso-offset" cube are offset by the same time for a given 

• 20. (dtnJ'To) ; pair^ ; ^^ y: 

By analysing the (dtn,; t 0 ) parameters in several 
picking times,' the method according to the invention can 
be used -in particular ' to' f determine the ' (V, r\) parameters 
necessary for processing including an anelliptic NMO 
25 correction of the traces of a CMP gather. 

The .said process includes the steps presented below 
in a simplified manner. 
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. During the preliminary step, an . analysis volume is 
defined including several nodes (dtn, t 0 ) . 

The following are performed, for all nodes in this 
volume: v 

■ r . - firstly,... according . to ^equation . .(,7a)., : the . static 
NMO correction of, traces in, the studied CMP. gather, as' a 
function of . the,, values ,, of . the ^.dtm,^ x 0 '. parang at the 
node considered, '" ^ static 'correction ' being valid 

for any picking time/ . V Y • . ... ■ 

- secondly, the semblance as a function of time 
associated with the correction made in the previous step 
is calculated. • • : 

Finally, the following are determined for each time 
in the several picking times: 

- firstly, the node (dtn, t 0 ) used to make an 
optimum correction, for example with regard to the 
semblance, criterion.; (this type of . criterion typically 
being ".' used ,' in. ■■ ief smicv^ • pidces^^ the 
horizontality" ' of : 'reflection' curves and to /determine- the 
'pliability of.' thei picking) ;/ I' , , .---J; > - ■ , 

- secondly, the/yalues' of the':, dtn ^' ,t 0 - parameters : at 
the said maximum semblance node v are convert: ed to: values 
of the velocity^? .^d';'/aneliip'ti6it^J-^' parameters at the 
said picking time considered:. - ^ . ■ . . ' , 

The velocity law (in . other words ' all (picking time , 
v >). pairs and the anellipticity law (all (picking time, 
T|) ) pairs are thus" defined. V : ^ 

Finally, the NMO correction for all seismic traces 
may be made using, these velocity V and anellipticity r\ 
laws in equation- (la) for the anelliptic offset 
hyperbola. 



Similarly, and as will be described in detail later, 
the (dtn/ t 0 ) parameters of the static PSTM migration 
(see equation 7b) ) may also be checked! The velocity V 
and anellipticity tj laws are then determined and may be 
used to make the PSTM" migration^ in' the^DSQR equation (lb) 
for the double' anelliptic offset hyperbola. V';"V 

Figure 2 illustrates the " effect "of the parameter To,, 
(and therefore the' anellipticity r\ according to equation 
(2) on the curvature residues after the NMO correction, 

. Note that the vertical, time' .. scale.' of the curve in 
Figure 2 is exaggerated so that this effect can be. 
clearly understood. 

Three curves are shown in Figure 2, for which the 
dtn parameter is fixed to the correct value and the 
parameter to is assigned different values. 

The central curve shows the case in which To. is 
eoual to its correct... value '.Soil: in other, words t when the 
corresponding equal toe its true, value 

horizontal. - _ \% . ^ s ; , * 1 : 1 V " ' * ? 

* - The upper curved shows , the. case in whi ch t 0 w i s.-' equal 
to .. a value- ■ t 02 :: less «;v than t .itsi . correct ... value ; Toi,: the 
corresponding' anellipticity " r\i> being greater than its 

true . value TJtrue • \ ' ';: ' y-X •' -A 

The lower' curve represents the v: case* .in which t 0 is 
equal to a value t 03 greater than its true value T 0 i, the 
anellipticity r| 3 being less than its true value T|true- 

Note that these lower and upper curves show that the 
"horizontality" of the corrected reflection curve is 
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acceptable at. near offset (x « 0) and at far offset (x » 

On the other hand, significant curvature residues 
are observed when, the offset x no longer tends towards 
5 one of these limiting values... 0 and . x««: In particular, 
■ particularly., significant residue^ f of -V an 

offset x centred" in ^ . 

For example, and as is shown, in Figure 2, when % 0 is 
equal to T02/ an RMOx 0 2 residual correction must be made 
10 to the of f set traces XrmoTo2- Similarly, when x 0 is equal 
to x 0 3/ a residual correction RMOt 0 3 must be made to 
offset traces x^o^- 

Due to these significant curvature residues, the 
offsets range can be used almost in its entirety in order 
15 to determine the anellipticity i\. 

Therefore, setting parameters in (dtn, . t 0 ) for the 
NMO correction make it possible to use available data for 
all;, of f sets:. ' (x^ between oZ and Xmax) ^ • in - determining^ the 
anellipticity T|. '* ' , * v \" \ x " , * ' 
20 As" has already been .ment i onedV • t hi s^is. • r no t - . the'- >: c as e • 

for NMO corrections' for which parameters' have been" set 
with V and rf for wfiich the estimate" of the an^eliipticity 
T| is made essentially using far offset data. v 
^ : Thus, v the" effect: of;^ 
25 is distributed on all of f set s^; " unlike the^ anellipticity r\ 
that. only affects far of f sets • Therefore/ the •"behaviour" 
of the parameter % 0 gives . a better constraint on 
anellipticity values. 



As already mentioned above, the _ optimum dtn and t 0 
parameters are determined within a 3D analysis volume 

(t 0 , dtn, x 0 ) . 

Several nodes (dtn, x Q ) , in other words several 
5 pairs, of" dtn, ..t 0 parameters f or . which values, are known, 
are considered in the, said analysis volume . ^ . 

Nodes are usually regularly spaced from each other,. 

by an increment Adtn on the dtn. axis and an increment At 0 
on the t 0 axis, 

10 Minimum values dtn^in, ^Tomin* to m in and maximum values 

dtn maX / Tomax/ toma* 1 of the' T 0 and t 0 respectively provide a 
means of defining the limits of the said analysis volume. 

Advantageously, plausible values of the velocity 
parameter V and the anellipticity parameter tj may be used 
15 to define a corridor .[dtn min (t 0 ) , dtn ma x(t 0 ) ] , [T 0 min(t 0 ), 
Tomax(tp)] inside, the said analysis volume.; •,; : .'■ 

This : corridor^:; restricts the v. "analysis , volume / and, 

therefore the number of "nodes (dtn, t 0 ) ' that have to be 

; r -\ : [ t .-\> ,., • • i . i ^ . % . •>-- ' . • •- ' -v }y y \ ^-t^^^s-'; yy^y 

considered to determine, the optimum pair (dtn, x 0 ) . .., v • 
20 • < ,;>••;, Although use of - the said corridor;. is.: 1 beneficial; for 
the ; efficiency of the process according to^ the invention, 

" it also provides a means of cohstrairiing v solution- towards 

; •'. :» * a 1 ; < y ; . \ . • 

the right phenomena, without .needing to : - consider 

incompatible (dtn, t 0 ) ,. (and therefore V, tj) ) pairs, for 
25 example related to multiple reflections or miscellaneous. 

interference phenomena V " -•-v • 

Figures 3a and 3b represent the said analysis volume 

(t 0 , dtn, t 0 ) in the context of an example of true seismic 

data processing during an NMO correction done according 
3 0 to the process according to the invention. 



Figure 3a represents three 2D panels a, b and c of 
the analysis volume: 

- panel a is a panel (dtn, x 0 ) with constant t 0 ; 

- panel b is a panel (dtn, t 0 ) with constant t 0 ; 
5 - panel c is a panel (x 0 / t 0 ) with constant dtn. 

Figure 3b diagrammatical ly shows'" the 3D: analysis 
volume (t 0 , dtn, x 0 ) and . three ' intersections with this 
volume along three planes with constant ' t 0/ to, s dtn 
i-espectively, each of these intersections being projected 

10 onto the corresponding side a, r b or c in Figure - 3a. 

Panels (dtn, t 0 ) with constant t 0 (panel a; in the 
above example) are panels in which the bispectral picking 
of the dtn and x 0 parameters is made, for example 
according to the maximum NMO correction semblance 

15 criterion, for the considered picking time t 0 . 

The tdtr^^ rto) , dtn max (t 0 ) ] \ [xomin (t 0 ) T 0ma x (t 0 ) ] 
corridor mentioned ■ above t for the effective ' analysis 
inside .. the analysis ^ volume • is also; shown in Figured 3a 
; . (lighter zone on each panel)' ..v^h: •. \"' V? -T': > \^. ^ A \ 

20 The seismic data sampling time step At 0 defines the" 

/■ difference between, tw^ 
for which the maximum associated semblance node (dtn/ x 0 ) 
is determined) and therefore the number of bispectral 
picking panels to be considered. 
25 Automatic picking:- also, provides a means . of 

extracting the dtn, t 0 parameters at a density that is 

greater when the At 0 increment between picking times is 
small . 

Sampling of the analysis parameters dtn and t 0 is 
30 directly related to the resolution of the . seismic 
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exploration; dtn and t 0 effectively have the same 
dimensions as seismic records. 

The systematic search for the maximum semblance, 
conventionally known in itself, provides a means of 
5 determining the (dtn,.. x 0 ) pair providing best focus, for 
a given picking time .to . . • . ..v ; ; •■ 

.... Parabolic interpolations about the .values, of nodes 
(dtn, T 0 ) can also provide a means of, evaluating values 
of the. dtn, x 0 . parameter between the different nodes 
10 that have actually . been/ picked./ And- in ; particular this 
type of , evaluation makes determination even ,more precise 
(in contrast to determination limited to nodes in the 
corridor) of the dtn, t 0 pai~ameter pair maximising the 
semblance function. 
15 Finally, the velocity V and - anellipticity r\ 

parameters are determined, always for the picking time t 0 
considered, by. using; the,, above/ mentioned, conversion 
equations (4) and - (5). : ... , •• :%( ■ .., 

;. Figure 4 shows the . comparison., between the. two 

20 approaches* (V, Van) and .(dtn, . x 0 ) , by representing the, 
:\ bi spectral picking panel for;r each,, for^a. given, picking 
-time...; : . .;; : '; { .' ... .; . ... ; 

The figure at the right illustrates the conventional 
approach (V, Van) for which the two axes are the velocity 
25 axes, (the anellipticity x\ being related to the ratio of 

■ " . . . . : . i(v a \ 

these two velocities according to r\ = ) . 

The left figure illustrates the approach (dtn, x 0 ) 
according to the invention for which the two axes are 
time axes. 
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It is important to note from the study in Figure 4 
that the dtn and x 0 parameters appear to be decorrelated. 
This "decorrelation" is striking when the two approaches 
^ . are . compared, spreading of. the (dtn,.. x 0 ) .. spectrum 
5 actually being much narrower than spreading of the, .(V, 
•Van) spectrum. . ,/v . . « ' '• . ; ., 

■ Therefore the picking* made in: the context . of - the 
(dtn, t 0 ) approach according . to the : invention . .is more 
precise than . the approach . carried . put in ' conventional 
10 methods..". . .: \.,/ v : , .;, r " . V >'- : ' 
Furthermore,. this : decorrelation enables filtering of 
the dtn and t 0 pickings separately,, while keeping NMO 
corrections. This is- not the case for the V and. t] 
parameters for which a reduction in one of the parameters 
15 must necessarily be compensated by an increase in the 
other and vice versa. , , ; 

.And due to interpolations and individual, filtering 

of parameters according to the invention, dtn and t 0/ it 
is then possible toy perform simultaneous interpolation 
20 and filtering, of the standard NMO correction parameters V 

; - ■ «■ and T] . ■* ""' r - • ■ * ■■ " - ^ T \ 

Figure 5 shows : iiori- linear -correspondence;;, according . 
' to equations (4) and-* (5) mentioned" above/ between the 
pair of time parameters,, (dtn, t 0 ), and the pair Vof 
25. velocity parameters- (V-V V^y. . ; - . ; v 

Figure. 6 shows, from left to right: 

- a CMP gather of real, seismic 'traces before- NMO 
correction; 

- the picking of the dtn parameter corresponding to 
30 this gather of traces; 
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- the picking of the x 0 parameter corresponding to 
this gather of traces. 

The straight line x 0 = t 0 on the picking of the 
parameter x 0 at the right in Figure* 6, corresponds to the 
5 purely hyperbolic reflection" curves • ' 

From right to.;; left, Figure , 7 represents the 
. semblance function :: : and- * the v velocity • V • and the 
1 anellipticity ' T] functions (see-' equation (4) and (5) ; ) 
deduced from ^ the : picking t of : " the dtn arid x 0 parameters 
10 shown in' Figure' 6. /. P^T-V.-v -V ' ' f' v * ; ■• :*■; 

On this example. pf v real seismic data processing/ it 
can be seen that the values of V and t| obtained generally 
corresponding to a semblance of more than 40%. 

The following description contains details of two 
15 particular embodiments of the invention. 

The first of these modes relates ; to a process for 
determining "optirnum parameters to make an ahelliptic NMO 
correction ; to the traces - of a CMP gather .; (see the 
diff exentV steps;" shown- inrth^ ^. y '--%; 

20 ' f»f*T With reference! to Figure. 8a^ this i first; embodiment 
includes an initialisation' step la during, which the 
, following operations are carried out in sequence: 

• determination/ of the.; limits .of the analysis- volume 

V-.. [ d tnminy ;. dtrimaxl'i 1 - [jbmin/ ^Omaxl [ t omia V > :: : t Omax ] x W^i •>.';'' 

25 ' ' •; calculation of /^O corrections CORRnmo . (equation 
(7a) ) for all offsets:- and,, for all, nodes . (dtn/ x 0 ) ' 
included in the analysis volume; 

- delimitation inside the analysis volume of the 

corridor [dtn min (t 0 ) , dtn max (t 0 )] [tomin (t 0 ) x 0ma x (t 0 ) ] of 
30 plausible velocity and anellipticity values. 
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Once the initialisation step la has been carried 
out, a step 2a to calculate velocity V(t 0 ) and 
anellipticity r|(t 0 ) laws is performed for each gather of 
CMP traces . 

5 * v ' This step 2a ' comprises": "r'rr.* ' ; 

- a' first "operation 1 ; 3 a : ' carried* out for each node 
(dtn, x 0 ) in the corridor defined in' -the \ .initialisation 
step ' la, during J vwhichr ; t are 
carried out' in sequence for- e¥ch picking time 1 to: : ^ r '' ;r ' 
10 ';■ - application of. static . NMO. corrections . CORRnmo . 

for all offsets along the 'corridor? precarculated during 
the initialisation step la; "> A - % 

- calculation of the semblance function on data 
corrected along the corridor using a time window 

15 appropriate to the dominant wavelet; 

■v -/summation "stack" ■ calculation) ; of 1 data 
corrected along the corridor (only near offset data can. 
be used advantageously, for, this purpose) ; i-^3 ,; '--'- 
: ■ ■ '\ ■ * .\-;> 4 a second, operation; t 4a carried, out for each pickings 
. 20 t ime . 1 0 ( the- said times being,, at intervals of , At 0 , between 
rtomin/ tomaxl during., which the following operations > are 
carried but : [ *i f:\ '£ ?U |? I i 

- search for maximum semblance in the [dt^m (t 0 ) , 
dtn ma x(to)] [Tomin (t 0 ) T 0 max (t 0 ) ] ; corridor of • the bispectral 

25 panel (dtn, t 0 ) ; ^ ■ " "T • ; , : r 

- check the fact that the position in "(dtn/ tqj "at* 
the maximum semblance corresponds to a summation extreme 
value for the same values dtn and t 0 ; 

- creation of the dtn(t 0 )i T 0 (t 0 ) and semblance 
30 ' (t 0 ) series; 
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- a third operation 5a designed to select and adjust 
the pickings obtained, during which the following 
operations are carried out: 

- increasing sort of the semblance series (t 0 ) ; 
5 - validation of pickings, dtn and t 0 for . which the 

distance in time * to f .the 'highest . semblance pickings is 
greater than. a. predefined value;., , . ; . 

; v . adjustment: ; of picking . dtn ^ and . x 0 }: : \ values ■ 

validated by., parabolic interpolations using surrounding 

10 .values;..;;,;;'-: ^;^",v".\.i^ :?i2 .y ' .i'.^i ': *;; v 4 & V '. "■ : *' 
",/?; - retention^, of- picked,:;, validated „ .and. adjusted 
values if it is possible to - calculate the Dix • interval 
.velocities; with . the., pickings .with, the highest 
similarities. . 

15 - a fourth operation 6a designed to use . equations 

(3). and (4) to convert picked, validated and adjusted and 
. ... retained . values of dtri'. and x 0 . during, operation;.. 5a, .=■■ into 
velocity V. and anellipticity r| " laws . / r " " : 

; The velocity V and" anellipticity y r\ " laws as a 
20 function ^ of time are' thus : perfectly determined . • The . 
. anelliptical.;-NMO ! " correction : of "seismic traces in the CMP 
gather can thus be done precisely. 

'.; The second" particular /; embodiment : of r the invention 
relates:: to:;" a'"tt; process;.,! f or determination' of ;, optimum . 
25 parameters- for anelliptical PSTM migration of traces in' a 
" • • CMP gather. :.. : ^v^':I^V V :.V :;: V > 'V ' 

This"-'.." second embodiment may ' be. included in a 
generalisation of the first mode discussed above. 

As has already been demonstrated, the use of the 
30 (dtn, x 0 ) parameters, enables static PSTM NMO corrections 
(see equation (7b) >. 
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In the context of a PSTM NMO correction, this 
application has the same advantages as previously 
discussed for the NMO obliqueness correction. 

More precisely, it will be noted that the first 
embodiment ..is. .5 only a. special " case': of " the second 
• embodiment • corresponding to the case of a zero migration 
aperture. ;•; . ' : //. "'• ' ' ** 

. /, With reference . to, Figure- 8b,' the second 'embodiment 
includes an- initialisation- step lb during which the 
following operations are carried out in sequence : 

' - determine the. i : limits of 1 " the /analysis volume 

- [dtn m i n , dtn^Mt] ',' Homin; i'omaxY, [tomin, t 0 ^] ; ? 

- calculate NMO corrections CORR PSTM (equation (7b)) 
for all nodes (dtn, x 0 ) included in the analysis volume 
and for all migration offsets inside the migration 
aperture; .: ; <v. . 

r delimitation -. : of . , ; plausible : f ': velocity ' ; and 
anellipticity values inside the' analysis volume of the 
corridor, [dtxw (to)"/ dtn ma *( t 0 ) ]« [x^ (t 0 ) ' W 1 (t* 0 ) K /' . 

Once this initialisation step lb': has been; done, the, 
said fir st embodiment t creates a. step. 2b to calculate the 
velocity V(t 0 ) "and anellipticity f| (t 0 Y, laws' for' each 
gather, of v . CMP traces . ., , §&4p «» 

This, step 2b comprises : * ..." ,///•:• /. ;.- . * ; 

- a first operation 3b done for each node (dtn, t 0 ) 
of the corridor defined "during the 'initialisation step 
la,., during which the following are / carried but in 
sequence : 

- for each offset class, operations for: 

- application on all midpoints inside the 
migration aperture along the corridor, of static 



corrections CORR PS tm precalculated during the 
initialisation step lb; 

- summation of corrected midpoints along the 

corridor. 

- for each. pi eking, time t 0 , operations to: 

- calculate, v; the : ■ semblance function' on 
corrected data along the corridor using a time window 
appropriate to the ^dominant;"; wavelet; .-Hir.- '■■ 

.. . - summation ./. (stack /-.. calculation):- of data 

corrected, along , the . corridor (only data with small 
offsets can be used advantageously for this purpose): 

- a second operation 4b carried out for each picking 
time t 0 , similar to the. operation 4a described. above, to 
create the dtn(t 0 ) , x 0 (t 0 ). and semblance t 0 series; 

- a third .. operation similar . to operation 5b 
described above, designed to select and adjust .the 
pickings obtained; ... . : ;• 

- a fourth operation; 6b similar : to operation- 5a 
described; above,;;; designed to., convert^ values^ of; dtn, x 0 
anto velocity Y and anellipticity ti laws. ; 

The - velocity ' , y r 'and , anellipticity rj • laws ' as **'a 
function of . time ..are thus- perfectly determined. And the 
PSTM migration of seismic traces in the CMP gather * may 
thus be done accurately . ..; . Lrv^fe'. ' : ;'>£,/•: :%l^v .^v - 

Figure 9 illustrates, the second embodiment" of the 
invention that has, just been described. V 

Seismic data are initially grouped into ■ iso-of f set 
cubes. 

As has just been described, NMO corrections CORR PSTM 
are applied (see first operation in step 3b) for each 
offset class (namely for each iso-offset cube) on all 
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midpoints. The arrow marked with the MIG label in Figure 
9 illustrates this operation. 

The midpoints' thus corrected are then summated 
during a second operation in step 3b, the arrow labelled 
5 STACK in Figure. 9, illustrating this operation; -V s l 

.... These two MIG and STACK operations are., specific; to 
the second, embodiment., of ...the invent ion J PSTM migration) . 

The following operations: : are l' performed '. in / the 
context of each of the two embodiments discussed .(the NMO 
' 10 ... obliqueness . . f correction, , as^/Jwe..^;. have,;,.; described, '. 
corresponding to the special case, of a PSTM migration 
with zero migration aperture). . : . . 

The semblance is then calculated for each picking 
time (second operation in step 3a for the NMO correction 
15 only, third operation in step 3b for the PSTM migration) , 
. the arrow labelled "semblance", in Figure. 9, illustrating 
this operation. ■ ;•■ * .•;'• 

. Obviously,' the MIG, STACK and semb lane e. calculation . 
operations, are implemented for each node (dtn,. t 0 ) . ..' 
. 20 The "automatic bispectrai : picking'! mentioned, in 

Fi 9ure ;; 9 corresponds to picking ' of maximum semblance 
parameters (dtn, t 0 ) for each picking time' t o: (operations ' 
4 a, 4b in Figures 8a, 8b respectively). The arrow 
labelled MAX illustrates , the search for the maximum 
25 semblance.' •• '" r "" ■= V ''^ ' 

■:■ ' : ; ^Finally, the pickings ' of 'the picked parameters' (dtn, 
x 0 ) are converted into velocity V and anellipticity r\ 
laws (operations 6a, 6b in Figures 8a, 8b and arrow 
labelled CONV in Figure 9) . 
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